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ABSTRACT 

• Weak lensing surveys are expected to provide direct measurements of the statistics 

of the projected dark matter distribution. Most analytical studies of weak lensing 
statistics have been limited to quasilinear scales as they relied on perturbative calcu- 
' lations. On the other hand, observational surveys are likely to probe angular scales 

less than 10 arcminutes, for which the relevant physical length scales are in the non- 
linear regime of gravitational clustering. We use the hierarchical ansatz to compute 
the multi-point statistics of the weak lensing convergence for these small smoothing 
angles. We predict the multi-point cumulants and cumulant correlators up to fourth 
order and compare our results with high resolution ray tracing simulations. Averaging 
over a large number of simulation realizations for four different cosmological mod- 
. els, we find close agreement with the analytical calculations. In combination with our 

work on the probability distribution function, these results provide accurate analytical 
Q\ [ models for the full range of weak lensing statistics. The models allow for a detailed 

■ exploration of cosmological parameter space and of the dependence on angular scale 

and the redshift distribution of source galaxies. We compute the dependence of the 
higher moments of the convergence on the parameters f2 and A and on the nature of 
gravitational clustering. 

o ■ 

' Key words: Cosmology: theory - large-scale structure of the Universe - Methods: 

' analytical 
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1 INTRODUCTION 



» I ■ Weak distortions in the images of high redshift galaxies due to gravitational lensing provide us with valuable information 
' about the mass distribution in the universe. The study of such distortions provides us a unique way to probe the statistical 
properties of the intervening large-scale structure. Traditionally, the study of gravitational clustering in the quasi-linear and 
non-linear regimes has been done by analyzing galaxy catalogs. However such studies can only provide us with information 
on how the galaxies are clustered. To infer the statistics of the underlying mass distribution from galaxy catalogs one needs 
a prescription for how the galaxies are biased relative to the dark matter. Weak lensing studies have the advantage that they 
can directly probe the statistics of the underlying mass distribution. 

Pioneering work in this direction was done by Blandford et al. (1991), Miralda-Escude (1991) and Kaiser (1992) based on 
the early work of Gunn (1967). Current progress in weak lensing can broadly be divided into two distinct categories. Whereas 
Villumsen (1996), Stebbins (1996), Bernardeau et al. (1997), Jain & Seljak (1997) and Kaiser (1998) have focussed on the 
linear and quasi-linear regime by assuming a large smoothing angle, several authors have developed a numerical technique to 
simulate weak lensing catalogs. Numerical simulations of weak lensing typically employ N-body simulations, through which 
ray tracing experiments are conducted (Schneider & Weiss 1988; Jarosszn'ski et al. 1990; Lee & Paczyn'ski 1990; Jarosszn'ski 
1991; Babul & Lee 1991; Bartelmann & Schneider 1991, Blandford et al. 1991). Building on the earlier work of Wambsganss 
et al. (1995, 1997, 1998) detailed numerical study of lensing was done by Wambsganss, Cen & Ostriker (1998). Other recent 
studies using ray tracing experiments have been conducted by Premadi, Martel & Matzner (1998), van Waerbeke, Bernardeau 
& Mellier (1998), Bartelmann et al (1998), Couchman, Barber & Thomas (1998) and White & Hu (1999). While a peturbative 
analysis can provide valuable information at large smoothing angles, it can not be used to study lensing at small angular 
scales as the whole perturbative series starts to diverge. 

A complete analysis of weak lensing statistics at small angular scales is not available at present, as we do not have 
a corresponding analysis for the underlying dark matter distribution. However there are several non-linear ansatze which 
predicts a tree hierarchy for the matter correlation functions and are thought to be successful to some degree in modeling 
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results from numerical simulations. Most of these ansatze assume a tree hierarchy for higher order correlation functions but 
differ in the way they assign weights to trees of the same order with different topologies (Balian & Schaeffer 1989, Bernardeau 
& Schaeffer 1992; Szapudi & Szalay 1993). The evolution of the two-point correlation functions in all such approximations is 
taken as given. Recent studies by several authors (Hamilton et al 1991; Peacock & Dodds 1994; Nityanada & Padmanabhan 
1994; Jain, Mo & White 1995; Peacock & Dodds 1996) have provided fitting formulae for the evolution of the two-point 
correlation function, which can be used in combination with the hierarchical ansatze to predict the clustering properties of 
the dark matter distribution in the universe. 

Cumulants and cumulant correlators are one and two point averages of higher order correlation functions. Scocciomarro 
& Frieman (1998) have shown that it is possible to predict the normalized non-linear one point cumumlants or Sn parameters 
in the highly non-linear regime by considering co-linear configurations of wave vectors in the perturbative expressions for the 
multi-spectra, a method labeled hyper-extended perturbation theory (HEPT). Munshi et al (1999c) combined this method 
with Bernardeau & Schaeffer's (1992) ansatz to predict amplitudes associated with distinct topological tree configurations of 
arbitrary order. They used this method to compute cumulant correlators in projected galaxy catalogs and found very good 
agreement. A similar analysis was done for the case of extended perturbation theory as proposed by Colombi et al. (1996) and 
was found to be in good agreement with that of HEPT. While such studies are generally done for projected catalogs it was 
shown by Hui (1999) that the hierarchal ansatze can also be used to predict the skewness of the convergence field associated 
with weak lensing surveys. 

In this paper we extend earlier studies of the lower order cumulants of the convergence field k{0q). We present a detailed 
theoretical analysis of cumulant correlators in the context of weak- lensing surveys (Munshi & Coles 1999). Using several 
realizations of weak lensing simulations for four different cosmological scenarios we compare the analytical results with 
numerical results for a range of smoothing angles and source redshifts. 

Several approximations are used in the analytical calculations of the higher order cumulants of the weak lensing conver- 
gence field k(0o). One such approximation is the Born approximation which neglects all higher order terms in the photon 
propagation equation. This means that the relation connecting the line of sight integration of the density field 8 with the 
convergence field k is only approximate. Perturbative calculations were used earlier to show that such correction terms are 
negligible for lower order cumulnats (Bernardea et al 1997). The effects of source clustering on the statistics of the convergence 
field have been estimated by Bernardeau (1997) using perturbative calculations. Such analyses are not possible in the highly 
non-linear regime as the entire perturbative series starts to diverge. The only way to check the approximations is to compare 
the results with ray tracing experiments. 

In section Section 2 we briefly describe the ray tracing simulations used to check the analytical predictions. Section 3 
presents most of the analytical results necessary for computing cumulants and cumulant correlators in the highly nonlinear 
regime. The comparison of the analytical and simulation results is made in Section 4. In section 5 we discuss our result in a 
general cosmological framework. The ray tracing experiments placed the source at unit redshift; in the Appendix we study 
how our results depend on the source redshift. 



2 THE GENERATION OF CONVERGENCE MAPS FROM N-BODY SIMULATIONS 

Convergence maps are generated by solving the geodesic equations for the propagation of light rays through N-body simulations 
of dark matter clustering. The simulations used for our study are adaptive P 3 M simulations with 256 3 particles and were 
carried out using codes kindly made available by the VIRGO consortium. These simulations can resolve structures larger than 
30h~ 1 kpc at z — accurately. These simulations were carried out using 128 or 256 processors on CRAY T3D machines at 
Edinburgh Parallel Computer Center and at the Garching Computer Center of the Max-Planck Society. These simulations 
were intended primarily for studies of the formation and clustering of galaxies (Kauffmann et al 1999a, 1999b; Diaferio et 
al 1999) but were made available by these authors and by the Virgo Consortium for this and earlier studies of gravitational 
lensing (Jain, Seljak & White 1999; Reblinski et al 1999). 

Ray tracing simulations were carried out by Jain et al. (1999) using a multiple lens-plane calculation which implements 
the discrete version of recursion relations for mapping the photon position and the Jacobian matrix (Schneider & Weiss 1988; 
Schneider, Ehler & Falco 1992). In a typical experiment 4 x 10 6 rays are used to trace the underlying mass distribution. The 
dark matter distribution between the source and the observer is projected onto 20 - 30 planes. The particle positions on each 
plane are interpolated onto a 2048 2 grid. On each plane the shear matrix is computed on this grid from the projected density 
by using Fourier space relations between the two. The photons are propagated starting from a rectangular grid on the first lens 
plane. The regular grid of photon position gets distorted along the line of sight. To ensure that all photons reach the observer, 
the ray tracing experiments are generally done backward in time from the observer to the source plane at redshift z — z s . The 
resolution of the convergence maps is limited by both the resolution scale associated with numerical simulations and also due 
to the finite resolution of the grid to propagate photons. The outcome of these simulations are shear and convergence maps 
on a two dimensional grid. Depending on the background cosmology the two dimensional box represents a few degree scale 
patch on the sky. Figure 1 shows a map of k for the LCDM model with a field size 3 degrees on a side and source galaxies 
taken to be at z = 1. For more details on the generation of /-c-maps, see Jain et al (1999). 

While resolution effects are important for small angular scales, finite volume corrections play an increasingly dominant 
role for larger angular scales. Although finite volume corrections have been extensively studied in the context of projected 
galaxy surveys (Szapudi & Colombi 1996; Munshi et al. 1999b) such results are not available for weak lensing surveys. The 
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Figure 1. Convergence or k map generated by ray-tracing through N-body simulations of the LCDM model. The field is 3 degrees on a 
side. 



angular size of our simulation box at z = 1 is about 3°. We therefore restrict our studies to smoothing angles less than 10' 
and use many realizations for a given cosmological model to estimate the scatter in the numerical results. These realization 
are constructed from subsets of the same simulation box by changing the direction of projection for the individual lens planes 
(for more details on the ray tracing simulations see Jain et al. 1999). 



3 CUMULANTS AND CUMULANT CORRELATORS 

In this section we provide the formalism for computing multi-point statistics using the hierarchical ansatz (Munshi & Coles 
1999). We use the following line element for the background geometry: 

dr 2 = ~cdt 2 + a 2 (t)(d X 2 + r 2 { X )d 2 Q.), (1) 

where the angular diameter distance is denoted by r(\) and the scale factor of the universe by a{t). r( X ) = K~ 1/2 sm(K- 1/2 X ) 
for positive curvature, r(x) = (— K)~ ll/2 sinh((— K)~ 1 ^ 2 x) for negative curvature and \ f° r zero curvature universe. For the 
present value of the Hubble constant, Ho, and of the mass density parameter, Qq, we have K = (f2o — l)#o- 



3.1 Formalism 

The statistics of the weak lensing convergence k is similar to that of the projected density filed. In what follows we will 
consider a small patch of the sky where we can use the plane parallel approximation or the small angle approximation to 
replace spherical harmonics by Fourier modes. The 3D density contrast along the line of sight when projected onto the 2D sky 
with the weight function lo(x) will provide us the projected density contrast or the weak-lensing convergence at a direction 7. 
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k(7i) = / d Xl uj(xi)S(r(xhi) (2) 



Assuming all the sources are at the same redshift (an approximation often used though not difficult to modify for a more 
realistic description), one can write the weight function as u)(x) = S/2ac~ 2 HQQ, m r(x)r(Xs — x)/ r (Xs)- Where Xs is the 
comoving radial distance to the source. Using a Fourier decomposition of 8 we can write 
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k(7i) = / dXi^(xi) / 7^—3 exp(ixifc|| + r8k ± ) <5 k , (3) 



(2^7 

where we have used ku and k± to denote components of the wave vector k parallel and perpendicular to the line of sight. In 
the small angle approximation, one assumes that k± is much larger than ku. 8 denotes the angle between the line of sight 
direction 7 and the wave vector k. 

Using the definitions we have introduced above, we can compute the projected two-point correlation function of k (Peebles 
1980, Kaiser 1992, Kaiser 1998): 

<K(7i)K(7a))c = r^^TT^ I 4\l °M8l) P(4t) , (4) 
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where we have introduced 1 = r(x)k±, a scaled wave vector projected on the sky. The average of the two-point correlation 
function smoothed over an angle 80 with a top-hat window W2(l0o) is given by, 
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<* 2 (7)> c =/ dXi^T / 7^ P(-77xW(W ). (5) 



Similar calculations for the volume average of the three-point correlation function and the four-point correlation function can 
be expressed in terms of integrals of the matter multi-spectrum B p . The volume averages for angular smoothing in effect 
smooth over a conical volume, which in the small angle approximation is cylindrical (Munshi & Coles 1999). The results for 
the three and four point function are: 

<^»« = f^/^ (6) 
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g / 'i h h h \ 
4 v r(x) ' r( X ) ' r(x) ' r(x) / J2 l i=o 

In our derivation of the above results we have assumed that both the smoothing angle #0 and the separation angle 812 are 
small. 

Cumulant Correlatos were introduced by Szapudi & Szalay (1997) as normalized two-point moments of multi-point 
correlation functions. Two-point cumulant correlators have already been measured in projected surveys such as APM by 
Szapudi & Szalay (1997). Using the measurements they were able to separate contributions from different tree topologies. 
The concept of two-point cumulant-correlators can be further generalized to multi-point cumulant correlators. 

We list below the results of our analysis for two-point cumulant correlators for the third and fourth order in n. The 
derivation of these results is very similar to their one-point counterpart; we have assumed that in addition to the smoothing 
angles being small, the separation angle between different patches is small too. This assumption is consistent with our 
assumption of the hierarchical nature of the correlation function in the highly nonlinear. The third order cumulant correlator 
is given by, 

<K 2 ( 7 i)M72))c = r ^r\d X [ d 2 h [ d 2 laTsf-lH'-rT'-rr)^ W2{h8 )W2{l28 )W2{h8 )^ P {il 2 8 1 2) (8) 
Jo r (X) J J Vr(x) r(x) r {X)^J2 l '=° 

In next section we will show that the normalized third order cumulant correlator C21 depends only on one hierarchal parameter 
Q3. Its fourth-order analogs C22 and C31 depend on two different hierarchal amplitudes, R a and Rb\ therefore the two fourth 
order cumulant correlators can be used to separate contribution from R a and Rt, while their one-point counterpart S4 can 
only measure an average contribution from the two different topologies. The two fourth order cumulant correlators are, 
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3.2 Hierarchical Ansatz 



In deriving the above expressions we have not used any specific form for the matter correlation hierarchy. The length scales 
that dominate the contribution for the small angles of interest are in the highly non-linear regime. Assuming a tree model for 
the matter correlation hierarchy in the highly non-linear regime, one can write for the most general case (Groth & Peebles 
1977, Fry & Peebles 1978, Davis & Peebles 1983, Bernardeau & Schaeffer 1992, Szapudi & Szalay 1993): 

— , — , {N ~ V) 

£n(ti, . . . risr) = 

Qn <<* H ^ ri ' r j)' ( U ) 

ck,N — trees labellings edges 

It is interesting to note that a similar hierarchy develops in the quasi-linear regime in the limit of vanishing variance 
(Bernardeau 1992), though the hierarchal amplitudes Qn,o are shape dependent functions in that regime. In the highly 
nonlinear regime there are some indications that these functions become independent of shape parameters, as has been proven 
by studies of the lowest order parameter Qz = Q (Sccociamarro et al. 1998). In Fourier space such an ansatz implies that all 
higher-order multi-spectra can be written as sums of products of the matter power-spectrum. 

Ba(ki,k 3 ,k8)^ Jk4=0 = 0(P(ki)P(ka) + P(ka)P(ka) + i J (ka)P(ki)) (12) 
fl 3 (ki,k 2 ,k3,k4)p fc . =0 = i? a P(ki)P(ki +k 2 )P(ki +k 2 +k 3 ) +cyc.perm. + i?bP(ki)P(k 2 )p(k 3 ) + cyc.perm. 



Different hierarchal models differ in the way they predict the amplitudes of different tree topologies. Bernardeau & Schaeffer 
(1992) considered the case where amplitudes are in general factorizable - at each order one has a new "star" amplitude and 
higher order "snake" and "hybrid" amplitudes can be constructed from lower order "star" amplitudes (see Munshi et al. 
1999a,b,c, and Melott & Coles 1999 for a detailed description). In models proposed by Szapudi & Szalay (1993) it is assumed 
that all hierarchal amplitudes of a given order are degenerate. 

We do not use any of these specific models for clustering and only assume the hierarchal nature of the correlation functions. 
Galaxy surveys have been used to study these different ansatze. Our main motivation here is to show that weak-lensing surveys 
can also provide valuable information in this direction, in addition to constraining the matter power-spectra and background 
geometry of the universe. We consider the third and fourth one point moments of k in terms of hierarchical amplitudes: 



<k 3 ( 7 )) c = (3Q 3 )C 3 [/4J (13) 
<k 4 (7))c = (12.Ro + m b )Ci[4 \ (14) 

(15) 

Equation (Jk|) was derived by Hui (1998). He showed that his result agrees well with numerical ray tracing experiments of 
Jain, Seljak and White (1998). The two-point cumulant correlators can be written as (Munshi & Coles 1999): 



(K,(7i)K«(7a))c = 2QzCz[Ke ng 12 \ + QlC 3 [i4 12 ] (16) 
(kJ(7i)k s (72)) c = &RaCi\K 2 o Ke 12 ] + 374C4[K£> Ke 12 ] + 6R a CA[^e ^e 12 \ + RbC4[nl 12 ] (17) 
(Ks(7i) k s(72))c = 4R b Ci[Kg Ke 12 ] + 47? a C 4 [«;e K9 12 ] + 4R b C4,[Ke o K0 12 ] + 4R a C 4 [>4 12 ] (18) 

(19) 

We have used the following notation to simplify our presentation and have neglected the smoothing corrections which were 
shown to be small by Boschan, Szapudi & Szalay (1994): 



K 0O = J d 2 \ P^-L.^W 2 {W ) (20) 
M 12 = J d 2 \ p(-L-^W 2 {W )exp(iW l 2) (21) 

CtK>e 12 ]=y r 2ftio( x ) K "o^i2 d X ( 22 ) 

The expression for Ct[Kg^Kg 12 ] can be computed numerically for specific models of the background cosmology and different 
values of the parameters m and n. F(x) denotes the various products of Ke a and K$ 12 which appear in the above expressions 
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Figure 2. Analytical predictions for the r.m.s variance of the convergence field k($o) are plotted against numerical ray tracing experiments 
as a function of the smoothing angle 6q. Solid, dotted, short-dashed and long-dashed lines correspond to theoretical predictions for the 
SCDM, LCDM, OCDM and TCDM models. Dots correspond to numerical computations averaged over several different realizations, 
with the error bars giving the scatter between different realizations. 

containing \ dependence. The values of the hierarchal amplitudes (which in general are insensitive to the background cos- 
mology) can be computed from numerical simulations or from hyper-extended perturbation theory (Scoccimarro et al. 1998, 
Scoccimarro & Frieman 1998). This will allow a detail comparison against simulation data. 

It is customary to define the normalized cumulants and cumulant correlator by the following equations while studying 
statistics of clustering of dark matter distribution: 
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In the context of weak-lensing surveys it is more convenient to define the normalized cumulants and cumulant correlators as 



Sn = 



<*"(7) )c 
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N/2 
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(K(7)2) (iV-l)/2 ( | Ks(7l)Ks(72) | ) l/2 

The Sn parameter defined above for weak-lensing studies are also denoted by Eat (Valageas 1999b) in the literature. 



(25) 
(26) 



3.3 Modeling the non-linear power spectrum 

In numerical computations of the variance and volume averages of higher order correlation functions, we have to include the 
effects of non-linearity in the evolution of the power spectrum. The evolution of density inhomogenieties introduces radial 
dependences in the expression for the power-spectrum. In linear theory it can be written as P(k, x) = a? (x)F 2 (x)P(k) = 
D\{x)P{k). Where D + {\) is the linear growth factor for the evolution of density perturbations and F(x) depends on the 
background geometry of the universe. Lahav et al. (1991) approximated it by the following fitting function: 



F(x) = ^ m a(x)" 1 [xf + hl^ix)- 1 + Q K 
x = 1 + Qiaix)- 1 - 1) + fi A (a 2 ( X ) - 1) ; / = ( 4V) 



(27) 
(28) 



8 




0.1 1 10 

d Q [in arc minutes] 

Figure 3. Analytical predictions for the normalized skewness 53 of the smoothed convergence field k(6o). Solid, dotted, short-dashed 
and long-dashed lines correspond to theoretical predictions for the SCDM, LCDM, OCDM and TCDM models. Dots correspond to 
numerical computations averaged over several different realizations, with the error bars giving the scatter between different realizations. 
A hierarchical form is assumed for the 3-point correlation function of underlying dark matter distribution. The amplitude Q3 is evaluated 
using hyper-extended perturbation theory (Scocciomarro & Frieman 1998). 
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Figure 4. Analytical predictions for the normalized kurtosis 54 of the smoothed convergence field k(8q). The lines and dots denote 
analytical and numerical results as in the preceding figures. The 4-point function is computed in analogy with the three-point function 
shown in figure 3. 
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Figure 5. Analytical predictions for the correlation function (k(7i)k(k2)) of the smoothed convergence field ft(6o) as a function of the 
separation angle 9±2 . The lines and dots denote analytical and numerical results as in the preceding figures. The smoothing angle is fixed 
at O = T. 
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Figure 6. The correlation function (k(7i)k(k2)) is shown as in figure 5, but with smoothing angle £>o = 2'. 



For o = l, ignoring the weak £7a dependence of the logarithmic growth factor, the above expression gives, 

F( X = 0) = ^n m (l + n° m 6 + ^ m -Q A )- 1 (29) 

On scales where the dimensionless power, A 2 (k, x)) = iTrk 3 P$(k, \), is comparable to unity nonlinear corrections become 
important. This is more important for models with low values of fl m . In the quasi-linear regime such corrections can be 
computed using perturbation theory. An alternative semi-analytical approach which maps the linear power spectrum to the 
non-linear spetrum was propsed by Hamilton et al (1991) and was subsequently developed by Nityananda & Padmanabhan 



1 III 


i i i i i i 

O /"^ 1 — \ ft A 

SCDM 




LCDM 


L -* *. 

\ 


OCDM 
TCDM 




■ 

\ — 
\« 

\ 

* \. 


I- 1- -it 


>^"» ^ 

i I\_ '-A 





= 1 * • 
■ 





0.1 1 10 

6 ]2 [in arc minutes] 



Figure 7. The third order cumulant correlator 021(812) of the smoothed convergence field k(8o) is plotted as a function of the separation 
angle 612. The lines and dots denote analytical and numerical results as in the preceding figures. The smoothing angle is 80 = 1'. 
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Figure 8. The third order cumulant correlator C2i(fi2) is shown as in figure 7, with smoothing angle 8q = 2'. 



(1991); Peackok & Dodds (1996), Jain, Mo & White (1995) and Padmanabhan et al. (1996)). This ansatz involves mapping 
the nonlinear spectrum for a given wave number k to the linear spectrum at a wave number &l which are related to each 
other by k L = k(l + A 2 (k, X ))~ 1/3 ■ The accuracy of the results depend on the function G which relates the nonlinear power 
on k with the linear power spectrum kL i.e. A 2 (k,x) = G[A L (kL,x)]- We use the fitting formula of Peackock & Dodds (1996) 
to describe the non- linear power spectra as a function of length scale and redshift. 
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Figure 9. The fourth order normalized cumulant correlator CiiiQxi) as a function of separation angle (?i2. The lines and dots denote 
analytical and numerical results as in the preceding figures. The amplitude Q4 for the four-point correlation function is evaluated 
using hyper-extended perturbation theory. We used Bernardeau & Schaeffer's (1992) model for the correlation hierarchy to separate the 
topological amplitudes R a and if,;,. 
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Figure 10. The fourth order normalized cumulant correlator 022(^12) as in figure 9, with smoothing angle 80 = 2'. 



3.4 Modeling hierarchical amplitudes in the non-linear regime 

As discussed above, different models for the hierarchical amplitudes predict different relations between the topological am- 
plitudes that contribute to a given order, e.g. for the fourth order correlation functions we have two different topologies, the 
"snakes" and "stars". Bernardeau & Schaeffer (1991) proposed that hybrid and snake topologies at a given order can be 
constructed from lower order topologies as the vertices of trees representing the correlation hierarchy obey a multiplicative 
relation. This will mean that R a — Q 2 . However note that the full prescription for Q, R a and other higher order "star" topolo- 



SCDM 
LCDM 
OCDM 
TCDM 




0.1 1 10 

12 [in arc minutes] 



Figure 11. The fourth order normalized cumulant correlator 031(612) of the smoothed convergence field k(8o) as a function of the 
separation angle 612- The lines and dots denote analytical and numerical results as in the preceding figures. The smoothing angle is 

e = 1'. 
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Figure 12. The fourth order normalized cumulant correlator 031(812) as in figure 11, with smoothing angle 80 = 2'. 



gies will need a more detailed modeling of the non- linear density distribution. However it was noted by Munshi et al (1999) 
that since at each order we have only one new star topologies and Sn parameters at each order are completely determined 
by the topological amplitudes of the same order, it is possible to determine these star topologies if we have a prescription 
to determine Sn parameters at various order. In the quasi-linear regime perturbation theory provides such a formalism and 
Bernardeau (1992, 1994) developed a scheme to predict the S„ parameters of arbitrary order in the quasi-linear regime using 
a tophat smoothing. Different extensions of perturbation theory have been suggested to predict the one-point statistics in the 
highly non-linear regime. Colombi et al (1996) proposed to replace the local spectral index n in the expressions for the Sn 



parameters with an effective spectral index n e //. However it was realized that while the same n e // can be used for all Sn 
for a given power spectrum, n e ff does not coincide with the local slope of the non-linear power spectrum. Sccocimarro and 
Frieman (1998) on the other hand noted that the values of the Sn parameters in the co-linear configuration in the quasi-linear 
regime describe the amplitudes for these quantities in the nonlinear regime in which they become shape independent. We 
use their hyper-extended perturbation theory (HEPT) to make concrete predictions for the Sn parameters for weak lensing 
statistics. Their model for the dark matter distribution predicts: 

Sa(n)=3 I i-yj r (30) 
- 54 - 27 2" + 2 3" + 6" 

S4{U) = 8 1 + 62^ + 33^ + 66" - (31) 

We find that our results for all Sn parameters at arcminute scales are very accurately described by the local spectral index 
n = —2.3 for all four cosmological scenarios. HEPT is very accurate in predicting one point cumulants; it is possible to 
combine it the Benardeau & Schaeffer (1992) ansatz to compute the topological weights associated with different type of 
diagrams as well (see Munshi et al. 1999 for more details). This is necessary to predict the cumulant correlators or other 
multi-point statistics. HEPT gives us Qz = 2.70 and Q4 = 9.52 for n — —2.3, and we get R a = 7.29 and Rb = 16.23 for 
snake and star topologies, respectively. As noted by Munshi et al (1999), one can also use EPT as proposed by Colombi et 
al. (1996) which gives us very similar results. Finite volume N-body catalogs introduce a bias in the determination of the Sn 
parameters (see.g. Colombi et al (1995)). In particular it has been noted that the Sn parameters are biased towards lower 
values if determined from finite volume catalogs. Elaborate schemes were developed to correct such finite volume effects and 
were tested extensively against numerical simulations (Colombi et al. 1995; Munshi et al. 1997). Such effects will also be 
important in the determination of various statistical quantities from ray tracing simulations, which use the outputs of N-body 
simulations. We have not included finite volume corrections in our analysis of lower order cumulants and cumulant correlators. 
A detailed study of such corrections will be presented elsewhere. 



4 COMPARING ANALYTICAL PREDICTIONS WITH NUMERICAL SIMULATIONS 
4.1 Cumulants 

Our study of one-pont statistics includes the variance, S3 and S4. We have computed these quantities as a function of 
smoothing angle ranging from 0.1' to 10'. We have used approximately ten realizations for each cosmological model studied 
by us, except for the LCDM model where only one realization was used. Sources galaxies were taken to be at z 3 = 1. The 
use of such a large number of realizations gives us the opportunity to probe the degree of fluctuation from one realization to 
another. We find that our results are very accurately described by analytical approximations. Larger scales are more likely to 
be affected by the finite volume corrections described above. Compared to two-point statistics, one point cumulants are less 
affected by the finite size of the catalogs. In our analysis of one-point statistics, we find very good agreement between theory 
and numerical simulations, which suggests that such effects are indeed negligible. It is also to be noted that we have used the 
small angle approximation in our analytical computations, which is expected to hold good for the small angular scales studied 
by us. 

In Figure 2 we plot the variance of the convergence field k 2 (Oo), smoothed with a tophat filter of smoothing angle 60 . The 
variance for the SCDM model is larger than in the other models. In most of the models the analytical predictions show a very 
good agreement with numerical results. Issues related to the effect of finite force resolution have been extensively dealt with 
in the study by Jain et al. (1999). We find that for all models a smoothing angle larger than 0.25' is sufficient to remove such 
effects. Compared to other models, the TCDM models shows a slight departure from the analytical results; this could be due 
to the presence of large-scale power in this spectra. In Figures 3 and 4 we have plotted the next order normalized moments, S3 
and £4. Again, the analytical models agree with the numerical results very well and it is also evident that the Sn parameters 
can directly probe cosmological parameters, especially Qo. At large smoothing angles one can see that analytical results start 
to overestimate the Sn parameters. This is related to the fact that we have considered the hierarchical amplitudes constant 
with scale; in practice they go through a transition from the non-linear regime to the quasi-linear regime. Nonlinear estimates 
of the Sn parameters based on EPT or HEPT therefore over-predict such amplitudes. Finite volume corrections also starts 
to become more important for larger smoothing scales. However despite such effects we find the agreement between analytical 
simulations and numerical results to be very good. A good match also indicates that the Born approximation is a very good 
approximation even for the smaller angular scales probed by us. 

Throughout our studies of cumulants we have used top-hat filters for smoothing the convergence field k. The analytical 
results we have presented here can easily be extended to other filters such as compensated filters, which may be more useful 
from an observational point of view. We hope to present results of such analyses elsewhere. 

Cumulants are normalized moments of the one-point probability distribution function (PDF) for the convergence field. 
Our analytical prediction maps cumulants of the underlying dark matter distribution to that of the convergence field. In a 
similar analysis it is possible to map the complete PDF of the density to the convergence PDF (Munshi & Jain 1999) 



4.2 Cumulant Correlators 



We have measured cumulant correlators from simulations for two different smoothing angles, 80 = 1' and 80 — 2'. For each 
do, several separation angles 812 were considered, ranging from from 812 = O.l' to 812 = 10'. We have measured C21 at third 
order and C31 and C22 at fourth order. For small separation angles when two patches overlap with each other, the cumulant 
correlators C pq become equal to the moments S p + g . The numerical results for C31 do correspond to S3 and C22 and C31 
correspond to S4 at small angular scales. 

In agreement with our findings for the one-point cumulants, we find that analytical approximations for cumulant correla- 
tors are very close to results obtained from numerical simulations. Figures 5 and 6 show the two point cumulant correlators, 
for smoothing angles 80 = 1', 2' respectively, while figures 7 and 8 show the third order correlators. The fourth order correlator 
C22 is shown in figures 9 and 10, with the same smoothing angles, and C31 is shown in figures 11 and 12. Cumulant correla- 
tors are more strongly affected by finite volume correction compared to their one-point counterparts. We find that for large 
separation angles #12, there are slight discrepancies between analytical and numerical results. The comparison can be more 
appropriately made by performing ray tracing experiments through bigger N-body boxes. It is also possible that for large 
angular scales, the corresponding physical scales are already in the quasilinear regime where the values of the hierarchical 
amplitudes are smaller than their values in the highly nonlinear regime. 

Since the cumulant correlators are normalized moments of the two-point joint probability distribution function, they are 
directly related to two-point statistics of the underlying fields. They can also be related to the biasing of dark matter halos. 
The results of such an analysis will be presented elsewhere. 



5 DISCUSSION 

Ongoing weak lensing surveys with wide field CCD are likely to produce shear maps on areas of order 10 square degrees. 
In the future, 10° x 10° areas are also feasible, e.g. from the MEGACAM camera on the Canada France Hawaii Telescope 
and the VLT-Survey-Telescope. Ongoing optical (SDSS; Stebbins et al. 1997) and radio (FIRST; Kaminkowski et al. 1997) 
surveys can also provide useful imaging data for weak lensing surveys. Such surveys will provide a very useful map of the 
projected density of the universe and thus will help us test different dark matter models and probe the background geometry 
of the universe. They will also provide a unique opportunity to test different ansatze for gravitational clustering in the highly 
non-linear regime in an unbiased way. Traditionally such studies have used galaxy surveys, with the disadvantage that galaxies 
are biased tracers of the underlying mass distribution. 

Most previous studies in weak lensing statistics used a perturbative formalism, which is applicable in the quasilinear 
regime and thus requires large smoothing angles. To reach the quasi-linear regime, survey regions must exceed areas of order 
10 square degrees. Since existing CCD cameras typically have diameters of 0.25° — 0.5°, the initial weak lensing surveys are 
likely to provide us statistical information on small smoothing angles, of order 10' and less. This makes the use of perturbative 
techniques a serious limitation, as the relevant physical length scales are in the highly nonlinear regime. We have employed a 
new technique based on the hierarchal ansatz to compute weak lensing statistics on these small angular scales. By combining 
the hierarchical ansatz with either hyper-extended or extended perturbation theory, we can make concrete predictions for one- 
point cumulants and many-point cumulant correlators (Munshi & Coles 1999). We found very good agreement between our 
non-linear analytical predictions and ray tracing results for angular scales ranging from sub arcminute to tens of arcminutes. 

Our studies also indicate that several approximations involved in weak lensing studies are valid even in the highly 
non- linear regime of observational interest. While a weakly clustered dark matter distribution is expected to produce small 
deflections of photon trajectories, it is not clear if the effect of highly over-dense regions capable of producing large deflections 
can also be modeled using a hierarchal ansatz and the weak lensing approximation. Our study shows that this is indeed the 
case; the effect of density inhomogeneities on very small angular scales can be predicted with high accuracy using analytical 
approximations. Our study confirms for example that the Born approximation is valid in the highly nonlinear regime. 

The two-point correlation functions and the projected power spectrum of the weak lensing convergence field k can be used 
to estimate the shape and normalization of the dark matter power spectrum. In this study we have shown that cosmological 
parameters can significantly affect higher order moments of the convergence field. Earlier studies have shown that the skewness 
of the convergence is sensitive to the cosmological matter density Q m in the quasilinear regime (Bernardeau, van Waerbeke 

6 Mellier 1997; Jain & Seljak 1997; Schneider et al 1998), as well as in the nonlinear regime (Jain, Seljak & White 1999; 
Hui 1999). We have extended the work on higher order moments to include multi-point moments and fourth order moments. 
The multi-point moments may be more easily measured from non-contiguous lensing surveys. We find that the qualitative 
dependence of the skewness of n on f2 m is present in these higher order statistics as well. In contrast, higher order moments 
of the underlying dark matter density fields are much less sensitive to background cosmological parameters (Jain, Seljak & 
White 1999). 

The finite size of weak lensing catalogs will play an important role in the determination of cosmological parameters. It is 
therefore of interest to incorporate such effects in future analysis. Noise due to the intrinsic ellipticities of lensed galaxies will 
also need to be modeled. The present study is based on a top-hat window function which is easier to incorporate in analytical 
computations. It is possible to extend our study to other statistical estimators such as M ap , which use a compensated filter 
to smooth the shear field and may be more suitable for observational studies (Reblinsky et al 1999; Schneider et al 1997). 
Our method of computing cumulants and cumulant correlators can also be generalized to study the bias associated with the 



convergence field. Such studies will provide an interesting method of measuring the statistics of collapsed dark objects (e.g. 
Jain & Van Waerbeke 1999). We hope to present such results elsewhere. 

It is interesting to consider an alternative to the hierarchical ansatz for the distribution of dark matter in the highly 
nonlinear regime. The dark matter can be modeled as belonging to halos with a mass function given by the Press-Schechter 
formalism, and spatial distribution modeled as in Mo & White (1996). When supplemented by a radial profile for the dark 
halos such a prescription can be used to compute the one-point cumulants of the convergence field. Reversing the argument, 
given the one point cumulants of the convergence field it is possible to estimate the statistics of dark halos. 

The one-point probability distribution function and its two-point correlations are the complements to the moment hi- 
erarchy studied in this work. Numerical studies of the one point probability distribution function (pdf) of the convergence 
field have been carried out by Jain et al.(1999) and a fitting function has been provided by Yang (1999). Recently Valageas 
(1999a) and Munshi & Jain (1999) have computed the pdf of k using the hierarchical ansatz and compared the results to the 
measurements from simulations. The success of their analytical results means that the pdf for a desired cosmological model 
can be computed as a function of smoothing angle and redshift distribution. Thus physical effects, such as the magnification 
distribution for Type la Supernovae, can be conveniently computed. Thus we have a complete analytical description, based 
on models for gravitational clustering, for the full set of statistics of interest for weak lensing - the one point pdf, two point 
correlations, and the hierarchy of higher order cumulants and cumulant correlators. This analytical description has powerful 
applications in making predictions for a variety of models and varying the smoothing angle and redshift distribution of source 
galaxies; numerical studies are far more limited in the parameter space that can be explored. 
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Figure 13. The variance (k 2 (0q)) as a function of smoothing angle 9q for two different redshifts. Note that the variance increases with 
source redshift whereas all Sn parameters decrease with redshift. 



6 APPENDIX 

In our ray tracing simulations we have placed the source galaxies at a redshift of unity. However it is interesting to check 
how these results change when source redshift is altered. We find that the variance and other one point cumulants are 
sensitively dependent on the source redshift. While the variance increases with source redshift for a given smoothing angle, 
the Sn parameters and the C pq parameters decrease with the increase in source redshift. For source redshifts z 3 < 1 the Sn 
parameters are very sensitive to change in source redshifts. At low source redshifts, lensing probes the dark matter in the 
highly non-linear regime whereas for larger redshifts successive layers of matter tend to make the probability distribution 
function of k more Gaussian, thereby reducing the values of Sn parameters. Nevertheless the ordering of the Sn parameters 
in different cosmologies does not change with source redshifts. For any given source redshift the Sn parameters are largest for 
open models and smallest for the TCDM model. Theoretical analysis by Valageas (1999b) has shown that for smaller source 
redshifts S3 oc z^ 1 ^ 2 and for larger source redshifts S3 oc (1 + Zs)" 1 . Our numerical results shows that this is indeed the case. 
Perturbative analyses of the dependence of the Sn parameters on source redshift give qualitatively similar dependences. 

In figure 13 we plot the variance for two different source redshifts z s — 0.5 and z 3 — 1. Figures 14 and 15 show S3 
and 54 for redshifts z 3 — 0.5 and z a = 2; comparison with the same quantities for z s — 1 plotted in figures 3 and 4 shows 
that the Sn parameters do decrease with source redshift for all smoothing angles. In figures 20 and 21 we have changed the 
source redshift continuously keeping the smoothing angles fixed at (9o = 1' (left panels) and 80 — 2' (right panels). The Sn 
parameters show a rapid variation for small redshifts whereas for larger redshifts they vary at a much slower rate. The two 
point analogs of the variance and Sn parameters i.e. the correlation function and C pq parameters shows a very similar pattern 
in their dependence on source redshifts. The two point correlation function ^(71)^(72) is plotted in figure 16 and the lower 
order cumulant correlators C21, C31 and C22 in figures 17, 18 and 19, respectively. For small separation angles, i.e. small 
values of #12, the cumulant correlators C pq match with their one-point counterparts, S p + q . 

In all of our plots we have assumed a specific form of the tree-level hierarchy as proposed by Bernardeau & Schaeffer 
(1992). The ansatz proposed by Szapudi & Szalay (1993) gives very similar results, as shown in figure 22. We have plotted 
C22 and C31 for z s = 1 (predictions for the third order cumulant correlator C21 are the same for these ansatze). 
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Figure 14. The normalized skewness S3 as a function of the smoothing angle £>o for two different redshifts. 53 is higher for lower source 
redshift due to smaller projection effects and greater nonlinear clustering of the lensing matter. 
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Figure 15. The normalized kurtosis parameter S4 as a function of smoothing angle 9o for two different source redshifts. Note that the 
ordering of the different curves does not change with source redshift. The Sm paraneters are largest for the OCDM model and smallest 
for the TCDM model for any given redshift. 




Figure 16. The correlation function (k(7i)k(72)) as a function of separation angle 812 for two different source redshifts z s . The smoothing 
angles are 80 = l',2' for the upper and lower curves at small 812, respectively. 
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Figure 17. The third order cumulant correlator C21 as a function of separation angle 812, for two different source redshifts z s . The 
dependence of the cumulant correlators C vq on redshift is very similar to their one-point counterpart. 
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Figure 18. The fourth order cumulant correlator C31 for two different source redshifts z s as a function of separation angle 8\2- The 
smoothing angles are 80 = l',2' for the upper and lower curves at small 812, respectively. The C pq parameters tend to increase for low 
source redshift like the Sn paramters. 
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Figure 19. The fourth order cumulant correlator C22 for two different source redshifts z s as a function of separation angle 612. 
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Figure 20. Variation of S3 with source redshift 2 S for fixed smoothing angle 60 = 1' in the left panel and £>o = 2' in the right panel. While 
the rate of increase is moderate for high redshifts the rate is much more rapid at smaller redshift. The growth of density inhomogeneities 
at lower redshifts causes the values of the Sjv parameters to increase, while for larger redshifts subsequent layers of lensing matter tend 
to make the probability distribution function more Gaussian. 
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Figure 21. Variation of S4 with source redshift z s for fixed smoothing angle 60 = 1' in the left panel and 60 = 2' in the right panel. 
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Figure 22. The analytical predictions for the cumulant correlators C22 and C31 are plotted using the Szapudi &; Szalay (1993) ansatz 
R a = Rh in combination with hyper-extended perturbation theory. The smoothing angles are 60 = l',2' for the upper and lower curves 
at small S12, respectively. We find that the predictions are almost identical to that of Bernardeau & Schaeffer's ansatz, although the two 
differ in the way they assign amplitudes to trees of different topologies of the same order. 



